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The electronic properties of shallow acceptors in p-doped GaAsjllO} are investigated with scan- 
ning tunneling microscopy at low temperature. Shallow acceptors are known to exhibit distinct 
triangular contrasts in Scanning tunneling microscopy images for certain bias voltages. Spatially 
resolved I(V)-spectroscopy is performed to identify their energetic origin and behavior. A crucial 
parameter - the STM tip's work function - is determined experimentally. The voltage dependent po- 
tential configuration and band bending situation is derived. Ways to validate the calculations with 
the experiment are discussed. Differential conductivity maps reveal that the triangular contrasts 
are only observed with a depletion layer present under the STM tip. The tunnel process leading to 
the anisotropic contrasts calls for electrons to tunnel through vacuum gap and a finite region in the 
semiconductor. 

PACS numbers: 71.55.Eq, 73.20.-r, 73.40.GK 
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I. INTRODUCTION 



III-V compound semiconductors have become a stan- 
dard host material for fundamental research. The ap- 
phcation of Scanning Tunneling Microscopy (STM) over- 
comes a major disadvantage of far-field techniques that 
inherently average over a large region of the sample com- 
pared to the lattice unit cell. Non-periodic local features 
on the smallest length scales can be probed^. Due to 




FIG. 1: (40x40nm^) STM constant current topography of the 
atomically flat GaAs (110) cleavage surface. The corrugation 
of the surface states is clearly visible. At the selected sam- 
ple bias +1.6V Zinc acceptors superimpose triangular shaped 
protrusions on the otherwise atomically flat surface. The tri- 
angles' heights vary depending on the depth of the dopant 
under the surface. The image shows acceptors buried up to 
eight monolayers. 



the GaAs's good cleavage properties^ even bulk defects 
and buried heterostructures can be studiecl^>^i^. Besides 
growing importance for future electronic devices (Accord- 
ing to the ITRS roadmap, gate lengths will reach the 
15nm range within five years^), this research field yields 
important information for modern solid state physics. 
Electronic properties of nanostructures close to the band 
edges can be studied with meV-resolution via Scanning 
Tunneling Spectroscopy (STS)ii^. 

In the past decade the electronic contrasts induced 
by shallow donors and acceptors in III-V semiconduc- 
tors have received extraordinary attention. FIG. [T] is 
a 3d visualization of a constant current topography ac- 
quired on a cleavage plane of a Zn-doped GaAs crystal. 
At certain sample bias voltages shallow acceptors show 
up as distinct triangular protrusions extending over sev- 
eral nanometers along the (100) directions of the crystal. 
The STM catches up triangles originating from accep- 
tors that are buried up to eight monolayers under the 
surface. Up to now anisotropic features are found for ev- 
ery investigated acceptor species (reports are known for 
2^9,10 Cdi^, Be^ and MnH). Additionally these 
features are observed in different III-V semiconductors, 
i.e. GaAsSiiiJ^ii^iii, InPi^ and InAsi^. Donors on the 
other side (e. g. Sii^iiiii^, Te^, and Sii2£) seem to show 
only circular symmetric features. They seem to comply 
with what is expected of an impurity hybridized with 
states of the nearly isotropic and parabolic conduction 
band. 

The origin of the acceptor's triangle with the dopant 
site located in the triangle's tip has been subject to sev- 
eral reports during the past years. At first it was sug- 
gested, that d-electrons could cause such an anisotropic 
shape or it could be the deformation of a circular contrast 
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due to tetragonal strain induced by the dopant atom* . 
Later it was argued that the triangle reflects an excited 
state of the acceptor wave functionii. 

Evidently topographic measurements do not suffice to 
resolve this puzzle. Therefore recent research is extended 
to spectroscopic measurements, i.e. spatially resolved 
I(V)-spectroscopy carried out on p-doped GaAs. I(V)- 
spectra on shallow acceptors show that the conductance 
with triangular shape is spread out over the whole band 
gap interval and continues to negative voltages, as well. 
Recently a description has been proposed for the occur- 
rence of this wide spread conductivity. If a depletion 
layer is present at the semiconductor surface, for small 
positive or negative bias voltages the only possible tun- 
neling channel involves tunneling through the depletion 
layer. A tunneling electron is exponentially decaying in 
the semiconductor, i.e., it is described by the complex 
wave solutions of the semiconductor band gap^^. For this 
transport channel the tunnel current is efficiently sup- 
pressed on the free GaAs surface because the depletion 
layer has a width of several nanometers. The acceptor's 
stationary negative charge locally perturbs the depletion 
layer and forms a double barrier potential. This results 
in a tunneling resonance, that enhances the tunnel cur- 
rent. The otherwise vanishing tunnel channel becomes 
visible. 

Within the proposed model the presence of a deple- 
tion layer is essential for the observation of anisotropic 
contrasts. This paper will focus on the dependence of 
tip induced space charge layers at the sample's surface 
with the applied bias voltage between tip and sample. 
The chosen system is Zn-doped GaAs. We discuss the 
differentiation of qualitatively different tunneling situa- 
tions. The results are based on the combination of mea- 
sured I(V)-characteristics and numerical calculations of 
the band-edge distribution under the STM-tip^*. En- 
ergy resolved conductivity maps (dl/dV-maps) allow to 
attribute the occurrence of the triangular acceptor con- 
trasts to specific tunneling conditions. In further anal- 
ysis every acceptor atom within the measured sample 
region is identified by integrated current maps. I(V)- 
characteristics of the undisturbed surface are discussed 
within this context. The TIBB(V) calculations and the 
dl/dV-maps point out that the observation of triangular 
features at the acceptors coincides with the presence of 
a depletion layer at the semiconductor surface. 



II. EXPERIMENT 

The experiments are performed in a custom built low 
temperature STM operating in UHV at a base pressure 
better than 2 x 10^^^ mbar. The system has a lateral drift 
of less than 0.5 A/h. This offers sufficient mechanical 
stability for high resolution I(V)-spectroscopy. The tips 
are electrochemically etched from polycrystalline tung- 
sten wire. Further processing is done in UHV by heat- 
ing, sputtering and characterization with field-emission 



before transferring the tip into the STM. Small adsor- 
bats that disturb the tunnel current are removed by 
millisecond bias voltage pulses during the measurement. 
This treatment results in reproducible highly stable tips 
which show no modifications of the apex atoms for sev- 
eral hours. The GaAs samples are cleaved in situ at room 
temperature and are transferred to the precooled micro- 
scope where they reach the equilibrium temperature of 8 
K within less than an hour after cleavage. The samples 
are conducting even at 4.2 K. The zinc doping concentra- 
tion of 5 X lO^^cm"^ establishes an impurity band with 
10 meV spectral width centered about 30 meV above 
the valence band edge--. At low-temperatures in bulk 
semiconductor, the Fermi energy is within this band and 
injected carriers are drained through it. 

The I(V)-spectroscopies presented in this work are ob- 
tained by stabilizing the tip at a certain setpoint tun- 
neling voltage and current and then switching off the 
feedback loop. The tip's position z(x,y) is recorded in 
the setpoint topography image. At the fixed tip sample 
distance z(x,y) first the apparent barrier height is mea- 
sured via dl/dz-modulation technique. The tip sample 
distance is periodically modulated by 0.5 A and the mod- 
ulation in the tunnel current is recorded. The apparent 
barrier height in eV is calculated under the assumption 
of exponential decay of the sample's density of States 
(DOS) into the vacuum and free electron mass of the 
tunneling particl o^^i^^ . Secondly an I(V)-characteristic is 
recorded by sweeping the sample bias voltage. The exact 
procedure of barrier height determination and numerical 
processing of the raw data will be treated in detail in the 
next section. 



III. NUMERICAL PROCESSING 

The measurement described above produces a data set 
which allows the extraction of spatially resolved apparent 
barrier height and energy resolved differential conductiv- 
ity maps of the area of interest. At first the topographic 
and spectroscopic data have to be disentangled: the tip 
sample distance has to be fixed at a certain tunneling set- 
point, the spectroscopic data is not recorded on a plane of 
constant height but on a corrugated surface. The height 
variations have impact on the I(V)-characteristics and 
would overlay energy dependent features^"*. This effect 
is minimized in the measurement by choosing a tunnel- 
ing setpoint where the topography is as smooth as pos- 
sible, but it cannot be eliminated, especially when I(V)- 
spectroscopy is performed with atomic resolution. With 
the knowledge of apparent barrier height and z-position 
of the tip, one can numerically project the measured I(V)- 
curves to a plane of constant height. To increase the 
signal to noise ratio a gauss-weighted gliding average al- 
gorithm is apphed to the I(x,y,V) data set. 

After the numerics the data set is displayed in different 
ways. I(V)-characteristics as seen in FIG.[2]give informa- 
tion of the overall voltage dependent behavior of certain 
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sample regions. They are obtained by averaging over all 
I(V)-curves in a specified range. The red curve in FIG. [2] 
for example corresponds to the red rectangle in the inset 
topography. The numerical derivative dl/dV of an I(V)- 
characteristic resembles the differential conductivity, i.e., 
the energy dependent variation of the sample's conduc- 
tivity (see lower part of FIG. [2]). dl/dV-maps show lat- 
eral variations in the differential conductivity at a se- 
lected sample bias. A common method to increase the 
signal-to-noise ratio in the band-gap region is to move 
the tip towards the sample during the I(V)-sweep. Usu- 
ally the changing transmission coefficient of the vacuum 
barrier is accounted for by a normalization of the dif- 
ferential conductivity dl/dV with a smoothed total con- 
ductance I /V^^ . In the following sections it will become 
evident that the I(V)-characteristics at the acceptors ex- 
hibit non-monotonic behavior. Therefore the following 
measurements are carried out without tip movement and 
dl/dV-maps are presented to avoid artifacts created by 
the normalization procedure. 



IV. LOCAL I(V)-SPECTROSCOPY 

Figure [5] presents the results of a STS-measurement 
performed on a zinc doped region. The upper part of 
the graph shows the I(V)-signal and as an inset the set- 
point topography. The lower part shows the numerically 
derived differential conductivity. At the chosen setpoint 
bias voltage 1.85V and tunnel current InA only faint cir- 
cular depressions of two buried zinc acceptors are visible. 
Thereby no artifacts originating from the height modula- 
tion above the acceptors occur in the I(V)-measurement. 
The red I(V)-characteristic is evaluated near one of these 
contrasts and marked by the red square in the topogra- 
phy. The blue curve represents an I(V)-curve evaluated 
above an undisturbed region of the sample marked by 
the blue rectangle. The data points represent raw data 
and the solid lines are smoothed to reduce the noise level 
of the dl/dV-signal as mentioned above. Both regions 
yield semiconductor characteristics with valence (nega- 
tive voltage) and conduction (voltage larger than -f 1.5V) 
band components separated by a band gap region of 
about 1.5V. Above the acceptor (red curve) two groups 
of peaks are visible with a spectral weight of about 10% 
of the setpoint current (lOOOpA). The tunnel current on- 
set is as low as 0.4V. This is most prominently seen in 
the dl/dV signal: While the red curve exhibits drastic 
modulations, the blue curve on the undisturbed surface 
shows nearly no conductivity over the whole voltage in- 
terval from OmV to 1550mV. Both conductivity peaks 
arc followed by negative differential conductivity (NDC). 
This points to the fact that the measurement is not only 
determined by the sample's local density of states. Usu- 
ally the tunnel current is assumed to be the integral 
over the sample's LDOS weighted with an energy de- 
pendent transmission coefficient which is monotonically 
dependent on the applied bias voltage^^. Therefore the 
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FIG. 2: Spatially resolved I(V)-spectroscopy. The spec- 
troscopy is measured above a part of the clean GaAs cleavage 
surface. The setpoint sample bias is -1-1. 85V and the tunnel 
current is set to InA. Two subsurface Zn acceptors are visible 
in the topography. They appear as faint circular depressions. 
No z-shift was applied. The red curve is evaluated in the 
small rectangle above the contrast of the lower acceptor and 
shows the energetic distribution of the acceptor induced con- 
ductivity. The blue curve is evaluated above the undisturbed 
region of the image (large rectangle). The open circles are 
raw data and the solid lines are the Gauss weighted average. 
The lower part shows the differential conductivity dl/dV of 
the upper two curves. Two intervals of negative differential 
conductivity are observed (NDC). The tunnel current above 
the acceptor rises up to 10% of the setpoint current. 



tunnel current should be a monotonic function of the 
sample bias. Within this approximation NDC is not ex- 
plained, so a more detailed model of the tunneling process 
is needed. 



V. LINKING ENERGY SCALES IN THE 
SAMPLE WITH THE APPLIED BIAS VOLTAGE 

Above all, the understanding of such measurements 
with 'anomalous' behavior requires exact knowledge of 
the potential distribution under the tip. The unpinned 
Fermi energy at the GaAs{110} surfaces'^^ causes the 
electric field between tip and sample to penetrate into 
the crystal. This induces bias voltage dependent space 
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FIG. 3: Numerically derived dependence of the tip induced 
band bending with the applied sample bias TIBB(V). The 
solid curve is evaluated for the measured value of 4.0eV tip 
work function whereas the dotted curve corresponds to 4.5eV 
work function. The main difference is the sample bias where 
the band bending vanishes (flat band condition). 
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FIG. 5: Histogram of the mean work functions of 25 STS- 
measurements on the sample system. The average value for 
all measurements is 3.96eV work function, assuming gaussian 
distribution of the data points. 
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FIG. 4: The tunneling barrier height is recorded for ev- 
ery pixel during the STS measurement via dl/dz-modulation 
technique. The map shows the lateral variations of the bar- 
rier height which correspond to the corrugation of the surface 
states. Histogramm analysis shows two peaks: 4.2eV above 
corrugation maxima and 3.8eV for corrugation minima. The 
weighted mean value is 4.05eV. 



charge regions under the tip. They distort the band edge 
ahgnment at the place where the measurement is per- 
formed. This effect is known as Tip Induced Band Bend- 
ing (TIBB)^^. 

The height and sign of the surface potential (TIBB) 
are evaluated by solving the Poisson's equation for the 
situation of a metallic tip in front of a semiconductor sur- 
face. For simplicity the calculations are performed one 
dimensional and follow the works of Feenstra et al.^^ and 
Koenraad et ali^. The results are plotted in FIG. [31 The 



following parameters are used in the calculations: The 
tip sample separation is approximated to 7A. This value 
is derived from previous works concerning the voltage 
dependent movement of the tip in z(V)-measurements^ 
and the z-movement of the tip, when a jump to contact 
occurs^i^^. As GaAs parameters the low temperature 
band gap of 1.52eV and electron afhnity 4.1eV are used^^. 
A p-doping concentration 5 x 10^^cm~^ and 31meV ion- 
ization energy for Zn are considered^"'^ . A crucial pa- 
rameter for the TIBB(V) dependency is the value of the 
tip's work function. Figure [3] shows two curves one eval- 
uated for 4.0eV the other for 4.5eV work function. The 
tip's work function directly enters the bias value of the 
flat band condition, so I(V)-spectroscopies are largely af- 
fected by modifications of the tip that change the work 
function. 

The actual value of the tip's work function for each 
spectroscopy can be derived from the measured appar- 
ent barrier height, if the electron affinity of the sam- 
ple is known. Because the sample surface is very flat 
and clean it is justified to assume that the local elec- 
tron affinity equals the {110} electron affinity of bare 
GaAs. Figure S] presents the spatially resolved appar- 
ent barrier height acquired simultanously with the I(V)- 
spectroscopy in FIGs. [2] & [SI The histogram of this bar- 
rier height map has two maxima at 3.8eV and 4.2eV re- 
spectively. They correspond to the work function val- 
ues on the maxima and minima of the measured surface 
states. The weighted mean value is 4.05eV. The tunnel 
barrier is approximated with a trapezoidal shape, like it 
is depicted in FIG. [Hid). Because the apparent barrier is 
measured by dl/dz- modulation, it is determined by the 
dominant transport channel at the setpoint topography. 
For this measurement (p-doped GaAs and -|-1.85eV sam- 
ple bias) it is electrons tunneling from the tip into the 
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conduction band of the sample. Electrons at {Ep + eV) 
have the highest transmission probability"^^, so they de- 
termine the measured barrier height. The trapezoidal 
barrier is related to an effective rectangular barrier using 
the Wentzel Kramers Brillouin (WKB) approximation. 
The measured barrier height of 4.05eV is reproduced 
for a tip work function of 4.0eV. It should be pointed 
out that the good matching of measured apparent bar- 
rier height and tip work function is only observed for 
the spectroscopies recorded on p-doped GaAs with set- 
point sample bias below about -I-2V, because then the 
dl/dz- measurement is done nearly at flat band condi- 
tion and the dominant tunneling channel is restricted to 
a narrow interval at the conduction band edge. When 
a spectroscopy is recorded at negative setpoint sample 
bias the measured barrier height exceeds 5eV although 
the tip work function is still 4.0eV. This is due to the 
enlarged tunnel barrier for electrons tunneling at the va- 
lence band. But still, if the band gap energy is included 
in the trapezoidal tunnel barrier, the tip work function 
can be extracted from these measurements, as well. 

A tip work function value of 4.0eV work function is 
typical for the tungsten tips prepared by the above men- 
tioned technique. Figure. [5] shows a histogram of the 
mean work function values of 25 STS measurements per- 
fomed on the p-GaAs sample system with several tips on 
different samples. The mean work function value of all 
STS-measurements is 3.96eV. It is worth noting that tip 
configurations with very low or very high work function 
values (e.g. 2.7eV or 4.7eV) show I(V)-curves with ab- 
normal behavior. Measuring the apparent work function 
is a good tool to validate the actual tip configuration and 
ensure reproducible imaging properties. 

In literature often the work function value of 4.5eV for 
bulk tungsten is used (see e.g. reff^i). But theoretical'^^ 
and experimentatfl studies reveal that this value strongly 
decreases when one moves from low index surfaces to 
higher index surfaces or increases the edge atom den- 
sity on the emitting W facets. In SEM measurements 
our tips exhibited typical apex radii of < lOnm. It is 
not surprising that the highly curved apex has a work 
function in the order 0.5eV lower than bulk tungsten. 
To emphasize the role of the work function for the in- 
terpretation of the STS-measurements the TIBB(V) in 
FIG. [3] is calculated for the measured value of A.OeV 
(solid curve) and for A.beV (dashed curve). Both curves 
demonstrate that the TIBB varies monotonically with 
the applied bias voltage and in each bias window. For 
depletion (for TIBB < OmeV) and accumulation (for 
TIBB >0 eV) respectively this dependence is approxi- 
mately linear. The slope of the curves depends on the 
charge density at the semiconductor surface that screens 
the tip's electric field'^'^. In the accumulation range the 
screening is dominated by free carriers and more efficient 
than in the depletion range where screening is accom- 
plished by the immobile acceptor cores. Therefore the 
TIBB(V) has a steeper slope for depletion compared to 
accumulation. The linearity reduces a bit if the energetic 



structure of the GaAs DOS is increased, but the situation 
does not change qualitatively. 

Although the overall shape of both curves is similar 
they differ strongly in the voltage position of flat band 
condition. For 4.5eV work function it is reached at 
1085mV and for 4.0eV at 1590mV. The point of flat band 
condition is crucial: This voltage point separates two 
qualitatively different tunneling regimes. At bias volt- 
ages higher than the flat band voltage the bands are bent 
upwards. A hole accumulation layer is formed under the 
tip and in the case of p-GaAs the highest states of the va- 
lence band become unoccupied. The valence band states 
become accessoble at positive bias. At smaller sample 
bias than the flat band voltage the bands are bent down- 
wards which results in a depletion layer. For electrons 
with lower kinetic energy than the band gap there are 
no unoccupied states at the sample surface. At positive 
voltage the only unoccupied states left are within the 
acceptor band in the bulk of the crystal. The only possi- 
ble tunneling channel requires that the electrons do not 
only tunnel through the vacuum barrier but also through 
the depletion layer in the GaAs crystal^^. Despite the 
change in the direction of tunneling for negative voltage 
the situation does not change until the sample bias over- 
comes the depletion layer potential. In the spatially re- 
solved I(V)-spectroscopy different characteristic features 
can be checked to validate the derived TIBB(V) and the 
flat band voltage position, i.e., onset of conduction band 
tunneling and observation of charge density oscillations. 



VI. SPATIALLY RESOLVED 
I(V)-SPECTROSCOPY 

The I(V)-spectroscopy presented in FIG. [5] was carried 
out on the whole sample region of(15 x 15)nm^. In FIG. [6] 
differential conductivity maps derived from this measure- 
ment are shown. Each map illustrates in grey scale the 
spatial variations of the differential conductivity at cer- 
tain voltages. In the presented region of (10 x 10)nm^ 
two subsurface zinc acceptors are visible. For reference 
the positions of the dopant atoms under the surface are 
indicated by white circles in each image. They are deter- 
mined by the center-of-mass of the circular topographic 
contraslji^. The topography is the flrst image of part (d) 
in the flgure (same topography as shown in FIG.[21). The 
dl/dV-maps are arranged with ascending order in bias 
voltage from bottom left to top right. 

The computed TIBB(V) relation (FIG. E]) shows that 
four different tunneling conditions are present in the mea- 
sured voltage interval. Accordingly, the dl/dV-maps are 
arranged in four columns FIG. [Hla), b), c) and d). The 
potential conflguration for a tunnel path perpendicular 
to the surface is depicted for each interval in a rigid- 
band model. The schemes show the alignment of the 
conduction band minimum (CBM), valence band maxi- 
mum (VBM) and the acceptor band (AB) with respect 
to the Fermi energy (Ej?). The bulk semiconductor is 
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(a) depletion, negative bias (b) depletion, positive bias (c) flat band, positive bias (d) acuumulation, positive bias 
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FIG. 6: Some representative dIdV-maps extracted from the STS-measurement shown in FIG. (2] The diflerential conductivity 
is shown in grey scale. Black is reduced and white enhanced conductivity with respect to the mean conductivity in each map. 
The black triangles in the maps for +1525mV to 1569mV are negative differential conductivity. The STS-measurement covers 
different tunneling situations. Each column of this figure has the same underlying potential configuration. They are sketched 
with rigid-band diagrams (description is given in the text): (a) depletion layer under the surface and electrons tunneling from 
the bulk of the sample through the depletion layer and the vacuum gap to the tip, (b) depletion layer and electrons tunneling 
from the tip through vacuum and depletion layer into the bulk sample, (c) flat band condition and electrons tunneling into the 
acceptor band and the conduction band, (d) accumulation of free holes under the surface and tunneling into the conduction 
band and the accumulation layer. 
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denoted (sc) and the space charge region (d). The right 
side of each scheme illustrates the approximate shape of 
the vacuum barrier (v) and the tip's Fermi energy plus 
applied sample bias voltage (eV) in the metallic tip (m). 
Filled states are grey and empty states white. Possible 
tunneling channels (It) are selected by the condition that 
we assume elastic tunneling only. They are indicated by 
arrows. The configurations are: 

(a) For negative sample bias voltages down to -135mV 
(FIGinia) the sample is in depletion. The surface poten- 
tial (TIBB) is negative and the band edges in the semi- 
conductor are pulled under the Fermi energy. Electrons 
are tunneling from the bulk of the semiconductor to the 
tip. The height of the depletion layer potential (TIBB) at 
the surface is larger than the applied bias voltage. At the 
surface no valence band states are present in the energy 
window of the sample DOS that can possibly contribute 
to the tunnel current (Ep to cV in the band edge dia- 
gram) . Electrons tunnel from within the sample through 
the vacuum barrier and also through the depletion layer 
in the semiconductor. 

(b) For positive sample voltages up to 1590mV (calcu- 
lated flat band condition TIBB(V=1590mV)=0meV, see 
FIG. 131) the depletion region persists. Despite the cur- 
rent direction is reversed, the situation stays the same. 
Until the sample bias overcomes the GaAs band gap 
(Egap = l-52eV) no unoccupied states of valence or con- 
duction band are accessible at the surface and electrons 
still have to tunnel through vacuum gap and depletion 
region (FIG. [Hb). Shortly before the flat band condi- 
tion is reached, electrons can access the conduction band 
{eV > Egap) and a second tunneling channel opens. 

The line (c) inicates the flat band condition at 1590mV 
in FIG. He). 

(d) For higher sample bias voltages the bands are bent 
upwards and the valence band maximum is pulled over 
the Fermi energy (FIG. (Hid). Due to the band bending 
free holes are accumulated under the tip. Electrons tun- 
nel directly into the conduction band and into the valence 
band. 

The dl/dV-maps are arranged according to this clas- 
sification. The overall evolution of the measured con- 
ductivity can be divided into two intervals: In the first 
interval triangular features are visible at the position of 
the dopant atoms. They are measured for the depletion 
regime at positive and negative sample bias (FIG. Ela) 
and b)). The triangles are positioned to the [001] side of 
the dopant atoms with the triangle's tip pointing to the 
[001] direction. In these triangular regions the conduc- 
tivity is greatly enhanced and subsequently followed by 
negative differential conductance as seen for the map at 
+1525mV. When the voltage approaches OV from either 
positive or negative bias the triangles fade slowly into 
the background noise but do not change qualitatively. 
At voltages next to OV (-19mV or up to -|-320mV) en- 
hanced conductivity seems to branch away from the lower 
acceptor atom. This is an artifact due to another accep- 
tor buried below the two prominent acceptors which will 
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FIG. 7: Current map derived from the dl/dV-maps represent- 
ing the amount of tunnel current originating from the sam- 
ple bias interval (-|-1547mV to -|-1668mV). The hole charge 
density oscillations are observed. The cross section directly 
through the center of the oscillatory pattern allows estimat- 
ing that the amount of tunnel current flowing into the hole 
accumulation layer is about rj — 27pA. 

be discussed in part (B) of this chapter. The second in- 
terval represents the observation of circular features at 
the acceptors. They appear at about -|-1580mV, which is 
at the computed flat band voltage position (FIG. [6]c)). 
Additionally the onset of strong surface corrugation is 
observed. For further increasing sample bias the circular 
feature becomes more prominent and shrinks towards the 
acceptor atoms' positions until it is no longer visible at 
+1756mV. 

A. Cross-Checking TIBB(V) with the 
Measurement 

The onset of the circular feature is related to an ac- 
cumulation layer under the tip and allows validate the 
calculated TIBB(V) dependence with the measurement. 
When flat band condition is passed, a few meV later the 
valence band edge is pulled over the Fermi energy. Free 
holes are accumulated under the tip. They are vertically 
confined by the extension of the accumulation layer. The 
lateral confinement is determined by the dimension of the 
tip (for our tips 10-20nm) and much weaker. Charge den- 
sity oscillations of the nearly free hole gas around charged 
defects are expected. Those oscillations have been pre- 
viously reported for electron accumulation layers on n- 
type semiconductor9i2i^i22ii£ and first indications of hole 
charge density oscillations (CDO) on p-GaAs have been 
observedi^. In the dl/dV-maps beginning at 1580mV a 
bright halo of increased conductivity surrounds the two 
acceptors. This halo begins with a diameter of about 
lOnm and decreases in size with increasing sample bias. 
This is the signature of a CDO in the hole accumula- 
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tion layeril. The TIBB increases when the sample bias 
is raised. The accumulation layer increases in energetic 
depth. More states fit into the tip induced quantum dot 
and the Friedel oscillation period decreases. 

It is possible that even when states exist in the sample 
they are not imaged because tunneling into them is sup- 
pressed by low drain rates into the bulk of the sample"^^ . 
The question whether the accumulation layer could al- 
ready be present at lower bias voltage should be ad- 
dressed. Because there is an additional tunnel channel 
into the conduction band for the voltage interval of in- 
terest it is worthwhile to deduce the fraction of the total 
tunnel current that is attributed to tunneling into the 
accumulation layer. A current map is derived from the 
I(V)-spectroscopy by integrating the differential conduc- 
tivity signal over the bias voltage interval from the onset 
of the CDO (1547mV) to the point where they are no 
longer visible in the dl/dV-map (1668mv). The inset of 
FIG. [7] presents this I(V)-map. The red line indicates 
the position of the plotted cross section. The corru- 
gation of the conduction band surface resonance over- 
lays the long-range signal of the CDO. The undisturbed 
surface (right part of the cross section) gives a signal 
of 157pA. The CDO causes an additional elevation of 
?7=27pA in the vicinity of the acceptor core. Thus the 
CDO has an impact of at least 17% on the tunnel cur- 
rent in the respective voltage interval. This estimation 
shows that a significant part of the tunnel current can 
be driven through the accumulation layer, proving that 
this channel is not suppressed and once a hole accumu- 
lation layer is formed, transport into it is well observable 
with the STM. The TIBB(V) calculation for 4.5eV work 
function would give rise to CDO already at about I.IV 
sample bias. For 4.0eV the CDO is expected slightly 
above -1-1. 6eV. Within the accuracy of the calculation, 
a measured onset of 1.58V fits well to 4.0eV tip work 
function. 

Another way to check the surface band alignment are 
the onset voltages of the different surface resonances on 
GaAs. According to the calculation for 4.0eV the onset 
of tunneling into the conduction band should begin at a 
bias voltage only a few meV lower than the GaAs band 
gap (about 1510mV) due to the negative band bending. 
This slight difference is not observable in the measure- 
ment with an effective resolution of about 20meV. If the 
tip's work function was 4.5eV the onset of tunneling into 
the conduction band should be observed at 1590mV. The 
onset of atomic corrugation with features of a conduc- 
tion band surface resonance indicates this onset in the 
dl/dV-maps. The energetically lowest surface resonance 
in GaAs on the conduction band side exhibits its cor- 
rugation along a {110} directioiti^i^. With our sample 
alignment it is imaged as white and dark lines running 
from top left of the dl/dV-maps to bottom right (best 
seen in the topography image). The onset is observed at 
1547mV, which is 43mV lower than expected for 4.5eV 
work function (see FIG. ??). This points to a work fimc- 
tion lower than 4.5eV. But the position and width of 



the surface resonances is not known to sufficient accu- 
racy. The onset of conduction band corrugation here 
only shows a tendency but is no proof. 

Concluding this section we find that the measured 
value of 4.0eV apparent work function is plausible. The 
TIBB(V)-curve calculated for 4.0eV tip work function 
is in agreement with characteristic features in the mea- 
surement. The comparison of the two TIBB(V) curves 
demonstrates that a variation of the work function largely 
influences the dependence of surface potentials with ap- 
plied bias, an experimental determination of this value 
improves the accuracy of the numerical calculations. 



B. Symmetry of the acceptor induced contrasts 

Now we focus on the voltage dependent features ob- 
served at the Zn acceptors in the dl/dV-measurement 
(FIG. [6]). The topography image shows no short range 
distortion of the atomic corrugation. Only two faint long 
range depressions superimposed on the atomic corruga- 
tion are visible. This indicates that the two zinc accep- 
tors are embedded under the two surface monolayers of 
the GaAs crystal. White circles indicate the positions of 
the Zn acceptor atoms. In the dl/dV-maps two qualita- 
tively different features are observed: (i) The aforemen- 
tioned CDO and associated circular features, (ii) Tri- 
angular features shifted along [001] to one side of the 
acceptor atom. 

In the maps from 1580mV to 1646mV the nearly circu- 
lar enhanced conductivity is visible and centered around 
the dopant atom's position. These feature have already 
been discussed in the above section and are attributed 
to charge density oscillations in the hole accumulation 
layer. For higher voltages only circular depressions re- 
main. They are the effect of the acceptor's stationary 
negative charge on the host crystal's density of states. 
Although the acceptor core is screened by free holes in 
the accumulation layer, residual charge remains for dis- 
tances smaller than the screening length. It is of the 
order of a few nanometers, here. The bands sketched for 
FIG. [6ld) are pushed further upwards by the negative 
charge. The conductivity locally decreases because the 
number of conduction band states available for tunneling 
decreases-. 

For small positive bias voltage triangular regions of 
enhanced conductivity appear at about 350mV and are 
most prominent on the conductivity peaks at 802mV and 
1387mV (compare to FIG. [5]). The enhanced conductiv- 
ity persists up to 1448mV and changes rapidly into nega- 
tive differential conductivity (NDC) as seen for -|-1525mV 
and -|-1547mV on the same triangular region. A similar 
NDC is found after the first conductivity peak, as well. 
For negative sample bias the triangular shaped enhance- 
ment becomes clearly visible at -58mV and persists un- 
til the end of the I(V)-measurement at -135mV. Both, 
enhanced and negative differential conductivity are re- 
stricted to a triangular region shifted to the [001] side 



of the accecptor atom with the triangle's tip pointing to 
the atom's position. This is true for the (110) and (110) 
surface. When the (110) or (110) surfaces are imaged, 
this triangle is reversed and appears on the [001] side of 
the acceptor^^'^^. The triangular region extends about 
2nm along [001] and the triangle's base is of comparable 
expansion. The contrast looks similar for both polarities 
(tunneling into or out of the semiconductor). 

The question whether the contrast symmetry does 
change when one passes from -58mV to +802mV is im- 
portant. Previously it has been suggested that the tri- 
angular contrast resembles the squared wave function of 
the first excited acceptor stateii. This approach is based 
on the assumption that for very small negative voltage 
the ground state, which is more elongated along [001] is 
imaged and that the contrast changes to the triangular 
one for lower bias voltage. The dl/dV-map for -19mV 
contains enhanced conductivity in the triangular region, 
but for the lower acceptor a faint branch of conductivity 
- hardly overcoming the noise level - seems to reach along 
[GDI] away from the acceptor atom. However, this branch 
is displaced by about Inm out of the symmetry axis of 
the triangular contrast. To get a conductivity image for 
very small positive voltage, the conductivity is averaged 
over the interval from OV to 320mV (see lower image in 
FIG.inib)). The resulting dl/dV-map looks similar to the 
one for -19mV. 

This apparent change in contrast symmetry for the 
"zero bias" conductivity maps is traces back to the over- 
lap of two acceptor contrasts. The prominent one orig- 
inates from an acceptor near the surface layers and the 
faint one from an acceptor buried deeper under the sur- 
face. Its contrast is also present in all other dl/dV-maps 
presented in FIG. [5] (e.g. for -58mV) but not always 
recognizable due to the images' color scale. An integral 
current map showing the tunnel current integrated from 
OV to 1.5V is used to identify all subsurface acceptors 
that influence the measurement. This map is shown in 
FIG. [51b). For comparison the topography is plotted, too 
(FIG. [51a). The sample area shown in FIG. [SI is marked 
with a dashed square. In the integrated current map even 
those acceptors are resolved that could not be identified 
in constant current topographies. As a guide-to-the-eye 
the identified acceptor contrasts are overlayed with grey 
triangles. The two acceptors nearest to the surface influ- 
ence the setpoint topography. They are visible as faint 
depressions. The integral current map resolves seven ac- 
ceptors in the region of the setpoint topography. Four of 
them are within the area of FIG. [HI 

According to the presented I(V) measurements there 
is no change in contrast symmetry observable for either 
small positive or negative voltage. In the depletion in- 
terval from -0.1 to +1.5V only triangular contrasts are 
observed near the acceptor atoms. In addition the direc- 
tion of tunnel current has no impact on these features. 
They persist over a voltage range of about 2V rather than 
being located at certain voltages. It is unlikely that the 
triangular contrast is solely attributed to the imaging of 
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FIG. 8: Integral current map of the STS in FIG. [2 The 
tunnel current is integrated over the voltage interval from OV 
to 1.5V and plotted color coded. The area inside the dashed 
rectangle in (a) is the same region as shown in the dl/dV-maps 
of FIG. [6] The integral current map allows the identification 
of all subsurface acceptors by their triangular contrast, even of 
those not visible in constant current topographies. The upper 
I(V)-curve in (c) is reproduced from FIG. [2l the lower one is 
an averaged curve above the undisturbed region marked in 
(b) with the white rectangle. The ordinate axis is magnified 
by a factor of ~20 with respect to the upper plot. 



a specific acceptor state's wave function, at least for the 
presented I(V)-spectroscopy. This is strong indication 
that additional tunneling processes are involved. 



C. I(V)-characteristics of the undisturbed surface 

The influence of deeply buried acceptor contrasts on 
the dl/dV-maps makes it interesting to examine the part 
of the sample where no separate acceptor contrast can 
be identified (white rectangle in FIG. [2]). The average 
I(V)-characteristic of this region is plotted in the lower 
graph of FIG. [He). The I(V)-curve recorded above the 
acceptor resonance is plotted in the upper part for refer- 
ence. In the band gap interval from OV to 1.5V there is 
faint tunnel current detected < hyA even on the undis- 
turbed surface. Comparison with the measured accep- 
tor induced conductivity reveals that the I(V)-curve pos- 
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sesses the same characteristic shape. We conclude that 
this is due to acceptors buried deep below the surface. No 
separate contrasts are visible, but rather a homogenous 
background of acceptor induced conductivity is detected. 
It is present on the whole surface, so the tunnel current 
is influenced by an acceptor at every point of the inves- 
tigated sample. This is consistent with the high doping 
density of the sample. 

D. Summary 

Bringing together the TIBB(V) calculations and the 
dl/dV-maps we find that the observation of triangular 
features at the acceptors coincides with the presence of a 
depletion layer at the semiconductor surface. The char- 
acteristic shape of the acceptor induced tunnel current is 
also present at the undisturbed surface, which indicates 
that this transport mechanism is active on every point 
of the sample. For bias windows when the anisotropic 
features are observed, the only possible tunneling chan- 
nel requires tunneling not only through the vacuum gap 
but also through this depletion layer. Neither conduc- 
tion band states nor valence band states are accessible 
directly at the surface. These observations link the tri- 
angular contrasts in the band gap window with a trans- 
port channel involving tunneling through a finite region 
in the semiconductor. When unoccupied sample states 
(bulk valence or conduction band or surface resonances) 
become accessible for tunneling at the surface, the trian- 
gular features disappear. This is measured as negative 
differential conductivity. The band related tunnel cur- 



rent shows only circular symmetric distortions at the ac- 
ceptors, that are caused by the acceptor core's negative 
charge. 



VII. CONCLUSION 

Spatially resolved I(V)-spectroscopy has been em- 
ployed to study the electronic properties of shallow accep- 
tors in p-doped GaAs. Against the background of several 
reports about the distinct triangular acceptor contrasts 
we focussed on the determination of the conditions in 
the sample that lead to these contrasts. Differential con- 
ductivity maps of two acceptors are analyzed with re- 
spect to the potential configuration under the STM tip. 
The triangular shape conductivity is spread over a large 
bias window of about 2V width. The band alignment for 
this bias range is derived on the basis of TIBB calcula- 
tions. It is the key finding that the sample is in deple- 
tion when the triangular shaped conductivity is observed. 
The tip's work function proves to be a crucial parameter 
for the TIBB calculation and in turn for the analysis of 
the I(V)-spectroscopies. The experimentally determined 
4.0eV work function for the tip apex is significantly lower 
than the often used 4.5eV for bulk tungsten. 
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